Oscillating solitons in a three component BEC 
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We investigate the properties of three component BEC systems with spin exchange interactions. 
We consider different coupling constants from those very special ones leading to exact solutions 
known in the literature. When two solitons collide, a spin component oscillation of the two emerging 
entities is observed. This behavior seems to be generic. A mathematical model is derived for the 
emerging solitons. It describes the new oscillatory phenomenon extremely well. Surprisingly, the 
model is in fact an exact solution to the initial equations. This comes as a bonus. 
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The idea of spinor condensates was first suggested in 
seminal papers of Ho and Ohmi [1]. The experimental 
creation of spinor condensates [2] , in which the spin de- 
gree of freedom, frozen in magnetic traps, comes into 
play, opened a new perspective to various phenomena 
that are not present in single component Bose Einstein 
condensates. These included the formation of spin do- 
mains [3] and spin textures [4]. A theoretical descrip- 
tion of the formation of spin domains can be found in 
[5]. A spinor condensate formed by atoms with spin F 
is described by a macroscopic wave function with 2F + 1 
components. Here we focus on the F = 1 case, which 
has been studied in a number of works [6, 7], Multicom- 
ponent vector solitons with F = 1 have been predicted; 
bright in [8] and dark in [9]. 



In this Letter we investigate the dynamics of an F = 
1 spinor Bose Einstein condensate for a wide range of 
scattering length. In particular we address the general 
physical problem of spin soliton collisions. For one spe- 
cific ratio of the scattering lengths, Wadati and cowork- 
ers in a series of papers [10] found a complete classifica- 
tion of the one soliton solution with respect to the spin 
states. They even presented an explicit formula of the 
two-soliton solution. One soliton solutions come in two 
classes: polar (spinless) solitons, and ferromagnetic soli- 
tons [10, 11]. Both can be generalized to a wider set 
of scattering lengths. In this Letter we focus on one of 
these classes - non-ferromagnetic (polar) solitons and an- 
alyze their collisions. As a result of these collisions we 
obtain new entities. These are solitons with populations 
oscillating between different spin components. Although 
the total spin is conserved, each emerging soliton carries 
nonzero spin. When we look closer we are able to describe 
these emerging entities by a model analytical expression. 



To begin with, we consider a dilute gas of trapped 
bosonic atoms with hyperfine spin F — 1. The wave- 
function in vector form is \& (x, t) = ^q, ^-i) T and 



it must satisfy the spinor Gross Pitaevski equation 
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Here the f a , (a = 1, 2, 3) are angular momentum opera- 
tors in a 3x3 representation and cq is negative to allow 
for bright soliton formation. We choose a natural set 
of scales characteristic for the problem: length scale - 

When we express 



and time scale t 
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Eq. (1) in these units, divide the equation by cq, all the 
coefficients but the ratio between self and cross nonlinear 
coupling which we denote by 7, will be equal to one. 

In the paper of Ieda et al [10] the authors considered 
this system with coupling constants C2 = cq = c < 
(7 = 1). In this case Eq. (1) describes a completely 
integrable system. The authors find N soliton solutions 
via the Hirota method [10]. In particular, they present 
both N = 1 and N = 2 solutions explicitly. 

Another special case occurs when 7 = and different 
spin components are separated. The two above cases are 
the only ones for which N soliton solutions have been 
found. To our surprise, when 7 differs from zero and 
one, as a result of colliding two N = 1 solitons a new 
class of stable, localized but oscillating entities emerged. 
The details of the numerical and analytical studies are 
presented below. 

Consider the polar one soliton solution to Eq. (1) 
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is parameterized by three angles r» , a» , fy. The ro- 
tation operator used here is lAi = U(j3i,on,6i) — 
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FIG. 1: Top: Density plot of 2|#±i| 2 (left) and \f \ 2 (right) 
for two soliton collision. Here 7 = — 1, ki = fca = 1/4, pi = 0, 
P2 — 1.65, ri = Ti. = /3i = /32 = 0, ai = and «2 = 7r/4. 
Bottom: Details of the initially stationary soliton pictured 
after the collision and to a different (much longer) time scale. 
Observe the oscillatory character of the soliton components 
and small momentum transfer (panels c and d). 

e -*/ 3 ft e -»/V e -*/ 3 ^ Note tnat q. i s ; n this case a 
dummy variable; our spinor (0,1, 0) T is not affected 
by e~ % ' 9i . Equation (2) describes a stationary mod- 
ulus soliton moving with velocity 2p. The norm (or 
total number of particles) of this soliton is: Nt = 
J dxn (x, t) — j dx . \j> = 4k, total momentum: 
P T = j dxp{x,t) = Jdx^ ■ {-id x ) * = pN T and the 
spin density: f a (x, t) = * t • f a ■ * = 0. This solution is 
valid for arbitrary 7. 

We collide two solitons of the type described by Eq. (2), 
with equal norms k, but otherwise different parameters. 
Due to the Galilean and rotational invariance we can set 
one of the solitons to be static with \ = (0, 1,0) T . The 
second soliton will move with momentum Ap and gen- 
eral on, Pi and Tj angles. In the case of 7 = 1, two 
colliding solitons of the kind described in Eq. 2 indeed 
emerge in the same form. Below we show that for dif- 
ferent configurations (values of 7) a new kind of solution 
can be observed after the collisions, and after an interval 
of radiation (Figs. 1 and 2). After the collision is com- 
pleted, we observe a pair of solitons that do not have a 
product structure. Spin and spatial degrees of freedom 
cannot be separated. The population oscillates between 
different spin components. We observe a well defined os- 
cillation frequency. We will refer to these new solutions 
as oscillatons. 

The result of the collision depends on the set of angles 
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FIG. 2: (Color online) Total energy of the soliton on the left 
in Fig.l as a function of time. Energy is lost to radiation. 
Energy in units of modulus of the energy of the initial soliton. 
Inset: peeks of modulus of |<&o| 2 (dashed blue line) and |\I/±i| 2 
(solid red line) of the emerging soliton. Observe the harmonic 
oscillations. 
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FIG. 3: Frequency of the oscillatons after the collision as a 
function of the collision parameter ai and relative momentum 
Ap. Note that there is a maximum for small momentum and 
a = 7r/4. The sector (tt/2,h) is a mirror image of this graph. 
Here 7 = — 1. 



(ai, /3j, Ti) and relative momentum Ap. In particular, we 
observe that the frequency of oscillations depends only 
on ai and Ap. The most striking feature of the collision 
is the appearance of a local spin density; the colliding 
partners carry away equal and oppositely polarized spin 
densities. Moreover the spin vectors will lie in the XY 
plane. The initial angle Pi determines the the final ori- 
entation, but not the frequency! This was verified in 
simulations. Figure 3 presents a three dimensional sur- 
face of the frequency of oscillations after the collision is 
completed, versus at and relative momentum Ap. The 
maximum frequency occurs at cc, = 7r/4. For a; = 
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we recover the 7 = 1 case; collisions are elastic. More 
interestingly, for aj = 7r/2 there is practically no interac- 
tion between colliding solitons. The oscillation frequency 
is also a function of relative momentum. We observe a 
clear maximum for small values of Ap. Also for small 
relative momentum we see in the numerical simulations 
that some energy is transferred to the external degrees 
of freedom. In our case the initially static soliton is set 
in slow motion (see Fig. lc and Id). 

We describe the collision in the reference frame of 
one of the initial solitons. The stability of the static 
wavepacket was investigated numerically. During the col- 
lision and for some time afterwards, part of its energy 
was lost due to radiation. Asymptotically, the energy 
of the stationary solution converged to a constant value, 
which was larger that the energy of the initial soliton (see 
Fig. (2)). 

In what follows we propose an analytic solution to 
Eq. (1) that fits the emerging oscillatons and turns out 
to be exact. We focus on the stationary soliton, since 
the moving partner is essentially described by the same 
wavefunction, but with a Galilean boost. We assume the 
following form of solution, hoping to describe the emerg- 
ing oscillatons 
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where Uf = U(f3f,af,6f). This ansatz, by its very form, 
describes oscillations between components. This may be 
seen by looking at the squares of the modulai 

|* 1 2 = 1 sin 2 a f (vl +f]- + 2f]+V- cos [cot -20/] )(4) 
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All the spin components oscillate with frequency u = 
/i+ — fi— [11]. The total global spin for this solution 
is a nonzero constant! Notice that a/ is related to the 
amplitude of the oscillations and by examining the spin 
density 



f = e a (*t . f° . *) = (r,l - rt) n, 
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where n = (sin a/ cos/3/, sina/ sin/3/, cosa/) T , we find 
that j3f influences the spin vector orientation. We see 
that Of only plays the role of a time delay. We also 
confirmed that the modulus of the local spin, |f| = 



V- 



— Tj^_ I , and the density 



arc time 



independent. 

Our anzatz substituted into Eq. (1) leads to the follow- 
ing system of two coupled ordinary differential equations 



-M± + (1 + 7)4 + (1 - 7)?? 



= 0. 
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The problem has thus been reduced to solving ordinary 
differential equations, a much simpler task than solving 
Eq. (1). Assuming that this can be done, our ansatz can 
be considered as an exact solution. 

We discovered that the results of the numerical studies 
reported above can be well described by our ansatz. For 
example, Fig. (4) compares predictions of our model with 
an emerging oscillaton found numerically. The match is 
perfect! The adjustment requires finding /j_ and the 
set of angles {/3/, a/, Of}, which is rather straightforward, 
and solving the coupled equations (7). 

The new solution (3) can be analyzed independently 
of the context of soliton collisions, as a new class of os- 
cillaton solutions to Eq. (1). A particular case is when 
7 = 1, for which a solution for the pair of equations (7) 
is obtained explicitly. For this 7 the equations for rj± are 
independent and each can be solved analytically. The 
solutions of interest arc 



r l± ( x ) = VM± sech (.y 
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To understand the significance of this result one should 
go back to Eq. (3) and appreciate its complex character. 
We found a solution that looks like being composed of 
two solitons! It oscillates with frequency lu = /i+ — 

For general 7, we have developed methods for dealing 
with two distinct regions. The first considers n + and 77 _ 
of the same magnitude, the second covers 77+ ^> r\- . 

When fi + = /_t_ we have no oscillations {uj = 0). This 
is equivalent to the solution given by Eq. (2). In this 
case 77+ = rj- = y/JI sech i^^/Jix) . When /i + and /_t_ are 
almost equal \x± = /z(l±5) we have found an approximate 
solution 



r}± (x) = a± seclV^ 1 (y/JIx) . 



(9) 



When we feed this ansatz into Eq. (7) we obtain two 
approximate equations for the amplitude 



/i± + Vmm± ~ (1 
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These equations yield the values of the amplitudes a± 
and thus conclude our search for an approximate solu- 
tion. We have checked this solution against numerics for 
values of 7 in the interval [—1/2, 1/2] and 5 between 0.05 
and 0.1. In all cases the fits were very close. 

On the other hand, consider the case when one of the 
etas is much larger than the other, say 77+ 3> r/-, any 7 
(excluding ±1). If we ignore n- in the first equation, we 
obtain for n + = y/2fi + /(l + 7) sech (^//ZZjTa;) . We next 
insert this value into the equation for rj- , where it plays 
the role of a potential: 

M- + 2 M+ (yTtO S ° Ch2 ^^"^ ^ 7? ~ = ° 

(11) 

Fortunately, this equation can be solved. The solution is 
?7_ oc sech £ (VmT.t) [tanh (y/jl^x)] , (12) 
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FIG. 4: (Color online) Comparison of the predictions of our 
analytical model with an emerging oscillaton found numeri- 
cally. Collision parameters are as in Fig. 1. We have adjusted 
H+, u) = fi + — fi- and the angles af,f3f,6f,Tf in Eq. (3) to fit 
the simulation. There is perfect agreement between the red 
(continuous) line corresponding to the analytic solution for 
|^±l| 2 given in Eq. (5) and the red dots obtained in the nu- 
merical simulation. Blue (dashed) lines and dots correspond 
to |* () | 2 (Eq. (4)). 



where e = y/ [i-/fi + , Pm" [£] is a Jacobi polynomial 
which must be finite of order n. This demands that //_ 
be one of the values given by 

^ = t(- 1 - 4ti+ V 1+8 it^) 9 ' (13) 

a spectrum. Our conclusion is that the values of //_ come 
in a discrete spectrum. 

In conclusion, we considered three component BEC 



systems with spin exchange interactions and a general 
coupling constant ci and any negative cq- We focused 
on two polar soliton inelastic collisions and discovered 
that robust oscillating entities are generated. Shortly 
after the collision, some of the energy is radiated away, 
but the emerging entities are again solitons, with popula- 
tions oscillating between spinor components with a well 
defined frequency and nonzero spin density. This fre- 
quency depends on the collision parameters. We derived 
a mathematical model which fits the emerging solitons 
perfectly. Surprisingly, the model is an exact solution to 
the initial Gross Pitaevskii equations. This solution is 
not limited to the simple product form of the spatially 
dependent and spinor parts. 

Experimental realization of the ideas presented here is 
perfectly feasible. It could be planned along the lines of 
the work concerning matter wave soliton trains [12], using 
Lithium. At first, a polar wavepacket should be produced 
in a tight optical trap, e.g. radial frequency 300 Hz and 
longitudinal frequency 100 Hz. Next, a magnetic field of 
700 G of constant magnitude should be applied to change 
the scattering length via Feshbach resonance to a slightly 
positive value, say 0.1 nm, and the longitudinal trapping 
frequency should be gently decreased to 10 Hz, to gener- 
ate an elongated structure. A pair of Bragg pulses can be 
applied to generate counterpropagating solitons. These 
solitons will travel towards the edges of the trap and be 
rebounded there. Before they collide, a specially pre- 
pared magnetic field, perpendicular to the original one, 
should tilt the relative spin direction. After the collision, 
which is inelastic, we measure the components of the os- 
cillatons. We numerically simulated this very sequence 
and confirmed the result. It was as expected. 

Bright solitons such as considered here exist only for 
attractive interactions. When interactions are repulsive, 
solitons, if they appear, are dark. They will be considered 
separately. 
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